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Abstract: The high energy neutrino detection by a km Neutrino Telescope placed in the Mediterranean 
sea provides a unique tool to both determine the diffuse astrophysical neutrino flux and the neutrino- 
nucleon cross section in the extreme kinematical region, which could unveil the presence of new physics. 
Here is performed a brief analysis of possible NEMO site performances. 



Introduction 

Neutrinos are one of the main components of the 
cosmic radiation in the high energy regime. Al- 
though their fluxes are uncertain and depend on the 
production mechanism, their detection can provide 
information on the sources and origin of the high 
energy cosmic rays. 

From the experimental point of view the detection 
perspectives are stimulated by the Neutrino Tele- 
scopes (NT's) constructed, like Baikal H] and 
AMANDA |2 |, or under construction like IceCube 
O under the ice and ANT ARES H in the deep wa- 
ter of the MediteiTanean sea. Here, also the exper- 
iments NESTOR [51 and NEMO |6| are in the R&D 
phase and, together with ANTARES, in the future 
could lead to the construction of a km"^ telescope 
as pursued by the KM3NeT project |7|. 

Although NT's were originally thought as de- 
tectors, their capability as Vr detectors has become 
a hot topic in view of the fact that flavor neu- 
trino oscillations lead to nearly equal astrophysi- 
cal fluxes for the three neutrino flavors. Despite 
the different behavior of the produced tau leptons 
with respect to muons in terms of energy loss and 
decay length, both and Vr event detection rates 



are sensitive to the matter distribution near the NT 
area. In principle, the elevation profile of the Earth 
surface around the detector may be relevant. In 
Ref. 1 8 1, some of the present authors calculated 
the aperture of the Pierre Auger Observatory tH 
for Earth-skimming UHE v^^, by using the Digi- 
tal Elevation Map (DEM) of the site (GTOPO30) 
[IJJ. In Ref. UJJ the DEM's of the under- water 
Earth surface, provided by the Global Relief Data 
survey (ET0P02) fT2l, was used to estimate the 
effective aperture for Vr and v^^ detection of a km"^ 
NT in the Mediterranean sea placed at any of the 
three locations proposed by the ANTARES, NEMO 
and NESTOR collaborations. In the present pa- 
per we further develop the approach of Ref. lITTI 
to evaluate the performances of a Mediterranean 
NT in the simultaneous determination of the neu- 
trino flux and the i^-Nucleus cross section in ex- 
treme kinematical regions (which may probe new 
physics, see e.g. [13]). Since the three different 
proposed sites for the under-water km^ telescope 
show event rate differences of the order of 20%, for 
the sake of brevity we report the results our analy- 
sis for the NEMO site only, which presents interme- 
diate performances. 
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Formalism and results 

Following the formalism developed in ||8] [TT| we 
define the km^ NT fiducial volume as that bounded 
by the six lateral surfaces So (the index a=D, U, S, 
N, W, and E labels each surface through its orien- 
tation: Down, Up, South, North, West, and East), 
and indicate with 17^ = {9a, 4'a) the generic direc- 
tion of a track entering the surface (see Figure 
4 of Ref. 1 1 1 1 for notations). We introduce all rele- 
vant quantities with reference to events, the case 
of i'^ being completely analogous. 

Let d<i>^/(di?^ dfia) be the differential flux of 
UHE i't + Dt- The number per unit time of r lep- 
tons emerging from the Earth surface and entering 
the NT through Sa with energy Er is given by 



dNr 

dt 



dna / dSa / d^;, 



d^,{E„na) 



dE„ dn„ 



J dE, cos {Bo) kl{E,,Er\ra,^a) ■ 



(1) 



The kernel kl{E^ , Er ;ra,^a) represents the pro- 
bability that an incoming crossing the Earth, 
with energy E^, and direction fta, produces a r- 
lepton which enters the NT fiducial volume th- 
rough the lateral surface dSa at the position fa with 
energy Er- For an isotropic flux and an exposure 
time T, the total number of t leptons (and similarly 
for muons) crossing the NT is 

Nr^T^jd^aj dSa j dE, j dEr 

Although the exact dependence of Eq. ^ on 
the neutrino flux and the neutrino-nucleon charged 
current cross section cr'^ may be quite compli- 
cated, basic physical considerations show that even 
a rough binning of the events for energy loss and 
arrival direction may be used to obtain information 
on both these quantities (see e.g. I14ill5l ). In par- 
ticular, in the following we shall consider the sum 
of the II and r contributions as the experimental 
observable, namely the energy deposited in the de- 
tector and not the energy and/or the nature of the 
charged lepton crossing the NT. In fact, only for a 
minor fraction of the detected events the nature of 
the charged lepton can be reliably established. 



According to Ref.s fT6l[T7l . the differential energy 
loss of the T leptons per unit of length in an un- 
derwater NT can be simply taken as dEr/d\ = 
-Pr ErQw, with = 0.71 X 10"^ cm2 g-1 and 
Pu, denoting the water density. Analogously, for 
muons one just needs to replace (3r with the corre- 
sponding value f3f^ = 0.58 x 10"^ cm-^ g^^. As- 
suming that the lepton energy loss in the NT by 
e.m. interactions, AEi, is just a small fraction of 
its energy at the entrance, Ei, we simply obtain 
AEi = \{fa,ila) Pi El Qw, where \{fa,^a) is the 
length crossed in the NT by the lepton whose track 
is defined the geometrical quantities fa,^a- 

Using these relations one can derive the spectrum 
of leptons detected in the NT as a function of their 
deposited energy, /S.E, and their arrival direction, 
= (0, 0), measured in the zenith-azimuth refer- 
ence frame 



d^iV 



d{AE)dn 



f^Sa [dE^ 

a=iJ.,T a 
1 d<i>,{E,) cos (0a) 



47r dE„ 



A(7-a,l^a)/3 



(3) 



By denoting with Xi a given bin in energy loss, 
and with Yj the one for the zenithal angle, we can 
integrate the expression (O to get the number of 
expected events in Xi x Yj, 



N,, 



E 

a.—fi.T a 



d{AE) / dn (4) 



dSa I dE.^^''^^'^'^^ 



47r dE„ 



\{fa,^a)Pa Qu 



To take into account the underwater surface profile 
one can numerically compute the above integral as 
described in Ref. JS) : by using the available DEM 
of the area near the NEMO site, one can isotrop- 
ically generate a large number of oriented tracks 
which cross the NEMO fiducial volume (see Figure 
4 of Ref. fTT]) and sample the above integrand. 
This technique allows also to account for the radial 
density profile of the Earth (we use the formula re- 
ported in fTSl). 

In order to study the sensitivity to both neutrino 
flux and cj'^ it is necessary to parameterize their 
standard expressions and the possible departures 
from them. In particular, we parameterize the flux 
as d$i,/di;^ df^a = C- 1.3- 10-8 {E^/GcVy^^ 
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Figure 1 : Angular distributions of + r) events 
collected in five years from a km'^ NT placed at the 
NEMO site (see text). 



Figure 2: {CT, D) region corresponding to the ob- 
servation of at least one event in each bin (for stan- 
dard cross section). 



GeV ^ cm ^ s ^ sr ^, which gives a standard 
Waxman-Bahcall flux QD for C = L» = 1. For 
the neutrino-nucleon cross section we use: 

{/ s 0.492 A 
(t) < E2 

, s 0.492 A , s 0.492 B 
(t) (t) E^ > E. 

where E\ = 10^'^ GeV is the energy below which 
the atmospheric flux is expected to dominate (so 
we consider only the region E > E\) and E2 — 
10^ ° GeV. In the low-energy bin this cross-section 
matches the standard expression ||20I for A — 1. 
A value of B significantly larger than 1 may be 
associated with new physics. Note that the factor C 
only enters via the product CT as a normalization 
and can be fixed to C ~ 1, considering instead the 
exposure time T as the effective variable. 

For illustrative purposes, in Figure[T]we report the 
event angular distribution, for a km'^ NT placed at 
the NEMO site in five years of operations. The solid 
and dashed lines correspond to events whose en- 
ergy loss in the detector belongs to the intervals 
lO^-^-lO^ GeV or > 10*^ GeV, respectively. The 
predictions are obtained for standard flux and cross 
section (^ = B = C = 1). In the plot are also re- 
ported the number of events Nij (see Eq. (HI) when 
we consider i = 1,2 for the previous two energy 



bins and j = 1,2 when the zenith arrival direction 
is between 0° and 90° or 90° and 180°. 

Clearly, for a very steep flux power-law index D, 
the number of events decreases. We shall re- 
quire that at least one event falls in each bin, in 
T years of running, in the case of standard cross- 
section; this rough criterion constrains the param- 
eter range that one experiment is able to explore to 
the brighter region of Fig. |2] corresponding to the 
intersection of the regions where Nij > 1, for all 
hj- 

As a preliminary result, in Fig. [3] we show the 
constraints (contours at the 68 % and 95 % CL) 
which can be obtained on the physical parameters 
A and D after the marginalization over C is made. 
Here we are assuming B = A, so that the plot 
represents the capability of the telescope to dis- 
entangle the energy dependence of the flux from 
the energy dependence of the cross-section (in the 
toy model where both are described by a single 
parameter). We performed a multi-Poisson like- 
lihood analysis f2Tl, in which the likelihood func- 
tion, L ~ exp(— x^/2), is defined using the fol- 
lowing expression for the i^ij being the event 
numbers of the reference model): 

= 2^ [(iV,, - N^^) + 7V° ln(iV° /TV,,)] • 

(5) 
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Figure 3: Marginalized contour levels in the 
{A, D) plane (for B) (see text for details). 

Conclusions 

We have performed an analysis of the capability of 
a km'^ NT in the Mediterranean to disentangle the 
high energy neutrino flux and neutrino-nucleon 
cross section in an unexplored kinematical region. 
Our statistical analysis exploits the dependence of 
observables on energy and arrival direction (under 
the hypothesis of an isotropic diffuse flux). Using 
a simplified toy model to parameterize fluxes 
and cross-sections, preliminary results confirm 
that this approach is very promising, and could 
potentially detect hints of new physics. Of course 
the real feasibility of such measurements will 
depend crucially on the size of the neutrino flux 
which fixes the time required to reach a reasonable 
statistics. A complete account of this research will 
be reported in a forthcoming publication. 
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